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ABSTRACT 

In this work, we investigate the rotational dynamics of the ginger-shaped near-Earth 
asteroid 4179 Toutatis, which was closely observed by Chang’e-2 at a distance of 
770 ± 120 meters from the asteroid’s surface during the outbound flyby (|Huang et all 
l2013a^ on 13 December 2012. A sequence of high-resolution images was acquired dur¬ 
ing the flyby mission. In combination with ground-based radar observations collected 
over the last two decades, we analyze these flyby images and determine the orientation 
of the asteroid at the flyby epoch. The 3-1-3 Euler angles of the conversion matrix 
from the J2000 ecliptic coordinate system to the body-fixed frame are evaluated to be 
—20.1°±1°, 27.6°±1° and 42.2°±1°, respectively. The least-squares method is utilized 
to determine the rotational parameters and spin state of Toutatis. The characteristics 
of the spin-state parameters and angular momentum variation s are extensively s tudie d 
using numerical simulations, which confirm those reported bv lTakahashi et al.l (|2013f) . 
The large amplitude of Toutatis’ precession is assumed to be responsible for its tum¬ 
bling attitude as observed from Earth. Toutatis’ angular momentum orientation is 

_l0 I r\ "1 cO 

determined to be described by Xh = 180.2 _q; 3 o and Ph = —54.75 _q; 3 oo, implying 
that it has remained nearly unchanged for two decades. Furthermore, using Fourier 
analysis to explore the change in the orientation of Toutatis’ axes, we reveal that the 
two rotational periods are 5.38 and 7.40 days, respectively, consistent with the results 
of the former investigation. Hence, our investigation provides a clear understanding of 
the state of the rotational dynamics of Toutatis. 

Key words: minor planets, asteroids: individual (Toutatis) - planets and satellites: 
dynamical evolution and stability - planets and satellites: interiors 


1 INTRODUCTION 

The Apollo-type near-Earth asteroid (4179) Toutatis was 
originally discovered on 10 February 1934 and remained a 
lost asteroid until it was once again detected by C. Pol- 
las and colleagues on 4 January 1989 in Caussols, France. 
From a dynamical viewpoint, the asteroid moves on an ap¬ 
proximately 4:1 resonant orbit at a large eccentricity with 
the Earth and has passed through a close encounter with 
the Earth every fou r years since 1992 dWhipple et al.iri993l : 
iKrivova et al.l[l994h . Dating back to the decadal years of 
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near-Earth flybys of the asteroid, the radar observations 
obtained from Arecibo and Goldstone reveal that Tou¬ 
tatis appears to b e an irregularly shaped asteroid with 


two distinct lobe s llOstro et ^ 1995: Hudson et al.l 1 19981 : 
lOstro et al.l 1 19991 . l2002l : iHudson et ~ 120031 ). Various types 
of ground-based observations indicate that Toutatis is a 
tumbling, non-principal axis (NPA) -rotating small body 
(iHudson et al.l[l9^ : IOstro et al.lll99Sll ). These effects, as ob¬ 
served in Earth-approaching flybys, have also been reported 


ments |Soencer et al.l 

19951: IHudson et al.ll 19951: lOstro et al.l 

1 19991: iTakahashi et al 

|2013|). 


The first near-Earth flyby for Toutatis occurred in De¬ 
cember 1992 at a distance of 0.242 AU, when the asteroid 
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again came into view. Optical observations were gathered 
from at least 25 sites around the world through an inter¬ 
national campaign. The observed rotational light curves of 
Toutatis appeared to be highly unusual, with a large ampli- 
tude and a non-perio dic long rotation period. Subsequently, 
ISuencer et al.l (|l995l l reported two major periods of complex 
rotation of approximately 7.3 and 3.1 days, as estimated 
from analysis of the data. In addition, they reported that 
Toutatis was the first asteroid to show strong photometric 
evidence of complex rotation. However, the authors did not 
clarify this complex rotation phenomenon. 

Furthermore, radar observations were performed by 
Goldstone in California and by Arecibo Observatory dur¬ 
ing Toutatis’ approach in 1992. The delay-Doppler images 
achieved a spatial resolution of 19 meters i n range and 
0.15 m illimeters per second in radial velocity. lOstro et al.l 
lll995h suggested a rotational period between 4 and 5 days 
b ased on these radar data. According to the investigations 
of iBurns et al.l lll97ll . [T973h . the damping timescale for the 
slow non-principal axis rota tion of Toutatis exc eeds the age 
of the solar system. Thus, lOstro et al.l lll995h noted that 
the spin state of Toutatis may be primordial. However, re¬ 
cent investigation suggests that YORP effects may slow the 
spin states of asteroids. Thus, Toutatis’ spin state remains 
a mystery. 

Based on these high- r esolut ion delay-Doppler radar ob¬ 
servations, [Hudson et al.l (Il995l f used a least-squares esti¬ 
mation to calculate Toutatis’ three-dimensional shape, spin 
states, and moment-of-inertia ratios. They showed that the 
dimensions along the three principal axes are 1.92, 2.40 
and 4.60 kilometers and that Toutatis rotates in a long- 
axis mode. The two major periods were found to be 5.41 
days for the rotation about the long principal axis and 7.35 
days for the long-axis precession about the angular momen¬ 
tum vector. The results derived from the ra dar data were 
incons istent with the solutions presented by ISoencer et al.l 
lll995ll . 

Moreover, [Hudson et al.l (Il 998ll ad opted the published 
optical light curves I Spencer et al.[[l 995ll and a rada r-derived 
shape and spin-state model ([Hudson et al.l [l995l f to esti¬ 
mate the Hapke parameters of Toutatis. The Hapke photo- 
m etric model wa s appli ed, and a X 2 minimization proposed 
by [Hudson et al.l l[l997lj was performed. The synthetic light 
curves that were generated based on their model provided a 
good fit to the optical data, with an rms residual of 0.12 mag. 
They showed that the combination of the optical data and 
radar observations led to an estimation of the spin-state pa¬ 
rameters for Toutatis that was superior to the radar-derived 
outcomes. The two parameters describing the moment-of- 
inertia ratios were determined to be 3.22 and 3.02, respec¬ 
tively. 

Bas ed on the t r iaxial ellipsoid shape and spin state 
given bv [Ostro et al.[ ([l995[ ~). [Krvszczvnska et all ([l999lf pre¬ 
sented the results of modeling the light curve variations of 
this unusual rotating asteroid by numerically integrating Eu¬ 
ler’s equation in combination with the explicit expression for 
an asteroid’s brightness as a function of Euler angles. They 
achieved good agreement between the observed and calcu¬ 
lated light curves. They emphasized that the light curves 
of Toutatis were dominated by the precession effect and by 
the superposition of precession and rotation, which resulted 
in an unapparent relationship between the rotation period 


alone and the light curves. This understanding yielded an 
appropriate explanation for the inconsistency between the 
rotational period of Toutatis determined from optical data 
([Spencer et_ah 1995[l and th at determined from radar obser¬ 
vations ( Hudson et al.lll995h . 

During the 1996 near-Earth approach, Toutatis was ob¬ 
served by the Goldstone 8510-MHz radar system. Based on 
the physic al model derived fro m the observations of the 1992 
approach, [Ostro et al.l ([l999l l analyzed the radar measure¬ 
ments and refined the estimates of the spin state of Toutatis. 
The combination of optical and radar data was proven to 
better predict the orientational sequence displayed in the 
images captured in 1996. After refinement, the two periods 
of Toutatis were updated and estimated to be 5.376 ± 0.001 
days for the rotation about the long principal axis and 
7.420±0.005 days for the uniform precession of the long prin¬ 
cipal axis about the angular momentum vector. These two 
parameters yielded moment-of-inertia ratios of 3.22 ± 0.01 
and 3.09 ±0.01. Thus, the orientation at the 2004 approach 
could be predicted in both inertial and geocentric coordinate 
systems. 

[Scheeres et al.[ ([2000h determined that mutual gravita¬ 
tional interactions between an asteroid and a planet or an¬ 
other asteroid can play a significant role in shaping the 
asteroid’s spin state. They analyzed the interactions of a 
sphere with an arbitrary mass and with Toutatis based on 
the radar-derived shape model. The results thus obtained 
could partially explain the phenomenon of Toutatis’ current 
unusual rotational state. It was demonstrated that the tum¬ 
bling spin state of Toutatis might have been caused by near- 
Earth flybys over its lifetime. This hypothesis enabled the 
estimation o f the mass distrib ution and moment-of-inertia 
for Toutatis ([Busch et al.[[201^ 1. thereby allowing the likely 
internal structure to be inferred. 

U sing radar observation s of five flybys from 1992 to 
2008, [Takahashi et al.l ([2013f l modeled the rotational dy¬ 
namics and estimated Toutatis’ spin-state parameters using 
the least-squares method. They calculated the Euler angles, 
angular velocities, and moment-of-inertia ratios as well as 
the center-of-mass (COM)-center-of-figure (COF) offset. By 
directly relating the COM-COF offset and the moment-of- 
inertia ratios to the spherical harmonic coefficients of the 
first- and second-degree gravity potential, they could de¬ 
termine the driving force of the external torque due to an 
external spherical body and evaluate the spin state. The 
terrestrial and solar tidal torques were considered in their 
dynamical models, and all aforementioned parameters were 
included in the variable state vector to be estimated in the 
study. Furthermore, the spin states and uncertainties were 
propagated to the 2012 flyby epoch. 

On 13 December 2012, the first space-borne close obser¬ 
vation of Toutatis was achieved by the second Chinese lunar 
probe, Chang’e-2, at a distance of 770 ± 120 {3a) meters 
from Toutatis’ surface ([Huang et al.l[2013a[ l. Optical images 
of the asteroid were acquired by one of the onboard engineer¬ 
ing cameras duri ng the outbou n d flyby . Through analysis of 
over 400 images, [Huang et al.[ ([2013a[ f estimated Toutatis’ 
osculating orbit, its dimensions along the major axes, and 
its orientations. The highest resolution of the images was 
better than 3 meters. New discoveries were made, includ¬ 
ing the presence of a giant depression at the large end, a 
sharply perpendicular silhouette near the neck region, and 
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direct evidence of boulders and regoliths. The geological fea¬ 
tures suggest that Toutatis may have a rubble-pile struc¬ 
ture. The physical length and width were determined to be 
4.75 X 1.95 km± 10%, respectively, and the direction of the 
+z axis was calculated to be (234.1°, 60.7°). They showed 
that the bifurcated configuration may indicate that Toutatis 
is of contact binary origin and that it is composed of two 
major lobes (head and body). 

In this work, we perform an extensive investigation 
of the optical images of Toutatis captured by Chang’e- 
2 , and we determine the orientation of the asteroid at 
the flyby epoch. In co mbination with radar observations 
dTakahashi et akl ll2013ll and references therein), we estimate 
the rotational parameters of Toutatis. Moreover, the solar 
and terrestrial tidal torques are considered in the establish¬ 
ment of the rotational dynamics model. The torque due to 
the misalignment of the center of mass and the origin of 
the b ody-fixed frame is evaluated to be insig:nifica nt at first 
order (iHudson et al.l[200l : lBusch et al.l[2'oi2l . l2014l '). Further¬ 
more, we incorporate the external gravitational tidal effects 
from the Moon and Jupiter in our dynamical model. 

Compared with the previous prediction 
dTakahashi et akl l2013l l. our results for Toutatis’ orien¬ 
tation, derived for Chang’e-2’s flyby epoch from both radar 
data and optical images, demonstrate good con sistency with 
the observational resul ts of the spacecraft dHuang et al.1 
l2013al : IZou et aP l2014l 'l. Our simulations reproduce the 
trajectory of the long axis in space, with a precession 
amplitude of approximately 60°. This high amplitude of 
Toutatis’ precession is supportive of its tumbling attitude 
as observed from Earth. The characteristics of the angular 
momentum variations is investigated in detail, and the vari¬ 
ation induced by the near-Earth flyby in 2004 is estimated 
to be 0.03%. The orientation of its angular momentum in 
space is found to be described by Ah = 180 . 2 ^q’ 3 o and 
/3h = —54.75j|iQ jQo, and therefore, this orientation has 
remained nearly constant over the past two decades. The 
rotational periods are estimated from the simulations to 
be 5.38 and 7.40 days for the rotation and precession, 
respecti vely. These values are in good agreement with the 
work of lOstro et al.1 dl999l l. 

This work is structured as follows: Section 2 presents 
the observational data, which comprise ground-based mea¬ 
surements and optical images acquired by Chang’e-2. In this 
section, we also analyze the optical data to derive the orien¬ 
tation of Toutatis at the flyby epoch. In Section 3, we model 
the rotational dynamics of Toutatis based on Euler’s equa¬ 
tion. The least-squares and multiple shooting methods are 
employed to fit the variable state vector and the correspond¬ 
ing results. The simulation results are presented in Section 
4. Finally, we conclude by discussing the innovations of our 
investigation compared with previous works. 


2 OBSERVATIONS 

As described above, Toutatis progrades on an approximately 
4:1 resonant eccentric trajectory with the Earth. Orbital 
determination and rotational parameters for Toutatis have 
been documented since the asteroid began to be continually 
observed in 1992. Since that time, ground-based observa¬ 
tions have been performed for its every near-miss of Earth. 


As is well known, on 13 Dec 2012, Chang’e-2 completed the 
first successful close flyby o f Toutatis and acqu ired numer¬ 
ous images of this asteroid (|Huang et al.ll2013al '). 

Using the released data from the Minor Planet Center 
and hundreds of optical observations from the ground-based 
observational campaign that lasted from July to Decem¬ 
ber of 2012, the orbital determination of Toutatis was pre¬ 
cisely achieved within uncertainties on the order of several 
kilometers, and the orbital parameters at the flyby epoch 
were calculated to be a=2.5336 AU, e=0.6301, i=0.4466°, 
D=124.3991°, u;=278.6910° and M=6.7634°. Hence, the ini¬ 
tial orbit can be integrated to calculate the relative positions 
of Toutatis with respect to the Sun, Earth, Moon and other 
major planets in the solar system, which are required for 
computing the external torques from the solar tides, the ter¬ 
restrial tides and the gravitational tides from other bodies. 

The positions of the major planets and the Moon are 
calculated based on the DE405 ephemerides released by JPL 
The gravitation of the Sun, the major planets and 67 
asteroids in the main belt as well as post-Newtonian effects 
are considered in the dynamical model to achieve the orbital 
integration of Toutatis. In addition, Chebyshev polynomial 
fitting is numerically implemented to obtain the position of 
the asteroid at any given epoch. 

2.1 Radar Measurements 

Radar measurements of Toutatis acquired by Goldstone and 
Arecibo from 1992 to 2 008 are used to so l ve for the asteroid’s 
rotational parameters. iTakahashi et al.l l|2013l ') presented 33 
sets of radar observations. Together with the orientation ob¬ 
tained by Chang’e-2 at the flyby epoch (see Section 2.2), 
we have 33 sets of ground-based observation outcomes, in¬ 
cluding Euler angles, angular velocities and one space-borne 
orientation parameter. The observational data are summa¬ 
rized in Table 1. The observational errors of the radar data 
are e stimated to be between 3° and 15° dTakahashi et akl 
l2013h for the Euler angles and between 2° day“^ and 10° 
day“^ for the components of angular velocity; these errors 
are taken into account in our fitting. 

2.2 Observations by Chang’e-2 

As mentioned previously, Chang’e-2 captured Toutatis’ sil¬ 
houette via one of the onboard engineering cameras at the 
time when the asteroid was approaching the Earth in De¬ 
cember 2012. The camera has a lens with a 54-mm fo¬ 
cal length and a 1024-by-1024-pixel CMOS detector. The 
field of view of the camera is 7.2° by 7.2°. The images 
were acquired during the outbound flyby of Chang’e-2 be¬ 
cause of the large Sun-Toutatis-Chang ’ e-2 phase angle o n 
the inbound route dHuang et al.l[2013al : IZhao et al.ll2014al i. 
The imaging of Tou tatis lasted approximately 25 minutes. 
iHuang et al.1 d2013al l reported the first panoramic image of 

^ ITakahashi et aO d201,'di used the DE430 planetary ephemeris 
which is more accurate for fitting the observation data of Tou¬ 
tatis. However, the position offset between DE405 and DE430 
is simply a few kilometers that will induce systematic errors of 
approximately 10 parts per million into our torque calculations, 
which is too small to change the conclusions of this work. 
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Table 1. Observational data for Toutatis from 1992 to 2012 ijTakahashi et aO l|2Q13h and references therein). 


Year 

Month 

Date 

Hour 

Minute Second 

Euler angles (°) 

Angular velocity (°/day) 

1992 

12 

2 

21 

40 

0 

(122.2, 86.5, 107.0) 

(-35.6, 7.2, -97.0) 

1992 

12 

2 

19 

30 

0 

(86.3, 81.8, 24.5) 

(-16.4, -29.1,-91.9) 

1992 

12 

4 

18 

10 

0 

(47.8, 60.7, 284) 

(29.1, -23.2, -97.8) 

1992 

12 

5 

18 

50 

0 

(14.6, 39.4, 207.1) 

(33.3, 8.2, -92.9) 

1992 

12 

6 

17 

30 

0 

(331.3, 23.7, 151.6) 

(6.6, 34.5, -95.8) 

1992 

12 

7 

17 

20 

0 

(222.5, 25.4, 143.9) 

(12.8, 25.4, -104.1) 

1992 

12 

8 

16 

40 

0 

(169.8, 45.5, 106.9) 

(-31.1 -21.9, -97.7) 

1992 

12 

9 

17 

50 

0 

(137.3, 71.3, 22.3) 

(11.8, -36.9, -94.9) 

1992 

12 

10 

17 

20 

0 

(103.1, 85.2, 292.6) 

(35.8, -8.9, -97.9) 

1992 

12 

11 

9 

40 

0 

(77, 85.7, 225.5) 

(31, 17, -96.3) 

1992 

12 

12 

9 

20 

0 

(42.8, 70.2, 133.2) 

(-1.3, 37, -95.9) 

1992 

12 

13 

8 

10 

0 

(13.7, 44.4, 51.9) 

(-38.3, 17.9, -97.3) 

1992 

12 

14 

7 

50 

0 

(323.7, 14, 0) 

(-70.5, -50.6, -91.1) 

1992 

12 

15 

7 

50 

0 

(193.2, 24.4, 21.4) 

(22.1, -26.6, -96.6) 

1992 

12 

16 

7 

10 

0 

(165.1, 46.4, 310.6) 

(33.4, -3.4, -93.7) 

1992 

12 

17 

6 

49 

0 

(130.6, 76.1, 234.9) 

(12.6, 33.9, -94) 

1992 

12 

18 

7 

9 

0 

(91.6, 81.6, 142.4) 

(-24.3, 29.6, -102) 

1996 

11 

25 

19 

48 

0 

(130.5, 78.9, 143.2) 

(-32, 16.4, -98.2) 

1996 

11 

26 

17 

51 

0 

(94.2, 88.1, 57.7) 

(-30.6, -18.7, -91.5) 

1996 

11 

27 

17 

34 

0 

(60.4, 81.2, 320.9) 

(10.7, -36.8, -94.7) 

1996 

11 

29 

15 

37 

0 

(349.3, 30, 168) 

(23.1, 28.9, -98.3) 

1996 

11 

30 

15 

47 

0 

(250.3, 14.2, 166.9) 

(-18.6, 32.1, -94.9) 

1996 

12 

1 

14 

23 

0 

(180.4, 37.6, 139.3) 

(-38.7, -0.5, -98.1) 

1996 

12 

2 

13 

43 

0 

(146.7, 64, 64.9) 

(-12.6, -34.8, -97.9) 

1996 

12 

3 

12 

20 

0 

(116.7, 81.4, 340.4) 

(24.3, -28.2, -98.1) 

2000 

11 

4 

17 

6 

0 

(110, 88.5, 30) 

(0, -32.5, -98.9) 

2000 

11 

5 

18 

1 

0 

(70.6, 84, 281) 

(34.5, -17.2, -97.9) 

2004 

10 

7 

13 

56 

0 

(79.19, 85.3, 365.2) 

(-2.5, -35.4, -109) 

2004 

10 

8 

14 

4 

0 

(44.9, 72.5, 263.1) 

(32.4, -18.1, -97.9) 

2004 

10 

9 

13 

57 

0 

(12.8, 47.3, 181.4) 

(29.7, 22.8, -98.1) 

2004 

10 

10 

13 

17 

0 

(327.7, 20.4, 124.1) 

(-10.7, 34.7, -97.3) 

2008 

11 

22 

10 

54 

0 

(119.5, 90.7, 92) 

(118.1, 90.4, 93.6) 

2008 

11 

23 

10 

45 

0 

(86.2, 85, 0.3) 

(-0.4, -36.2, -98.9) 

2012 

12 

13 

8 

30 

0 

(-20.1, 27.6, 42.2) 




Figure 1. Images of Toutatis. a: First panoramic image captured by Chang’e-2. b: Best-matching attitude of the radar model with the 
optical image shown in Fig. la. c: Illustration of the graphical frame of the camera, d: Combination and comparison of the radar model 
with the optical results. 
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Table 2. Comparison of optical and radar image results. 


Property 

Optical Image 

Radar Image 

Slope 

1.386 

1.340 

Ratio of length to width 

2.52 

2.08 


the asteroid in the sequence, which was acquired at a dis¬ 
tance of 67.7 km from Toutatis at a resolution of 8.30 m, as 
shown in Figure la. 

2.2.1 Attitude matching 

The 3D shape model derived from delay-Doppler radar 
imaging, as presented in Figure lb, is used to discern the 
attitude of Toutati s. The latest radar-derived shape model 
llBusch et al.l l2012ll was constructed based on additional 
radar measurements performed by Goldstone in 2000 and 
by Arecibo in 2004 and 2008, and in this work, this model 
is employed in combination with the optical images from 
Chang’e-2’s flyby to match the spin state of Toutatis. 

In general, the attitude of a rigid body in space can 
be determined from its rotations about the three axes of an 
orthogonal coordinate system. As Figure 2 shows, to obtain 
the attitude of Toutatis, we suppose that the axes /i, I 2 and 
I 3 are defined as follows: the mutually perpendicular axes 
l\ and I 2 extend through the center of Toutatis’ shape and 
along the directions of the long short axes in the image, 
respectively. In addition, I 3 is perpendicular to the image 
plane through the intersection of h and I 2 , and thus, these 
axes form a right-handed coordinate system. 

Considering t he render and the orientation of the cam¬ 
era’s optical axis dZhao et al.ll20143 l. the 3D radar-derived 
shape model of Toutatis can be rotated at an interval of 
1° for each of the three Euler angles about the three prin¬ 
cipal axes of its body-fixed frame to match its attitude to 
that shown in the optical images. As shown in Figure 2 , we 
choose three criteria related to the optical image to deter¬ 
mine whether the rotated model is consistent with the opti¬ 
cal results from Chang’e-2’s view direction: (1) the slope of 
the long axis, represented by the red line in Figure 2a; (2) 
the ratio of the long axis to the short axis, indicated by the 
green line; and (3) the obvious topography on the neck area 
connecting the two major lobes of the asteroid, as shown in 
Figure 2b and 2c. These features of the optical image can 
be reproduced by rotating the 3D radar shape model about 
its three principle axes. Figure lb shows the best approxi¬ 
mation of the attitude of the radar model to that indicated 
by the optical image shown in Figure la. The two images 
are quantitatively compared in Table 2. 

Between the optical image and the model, the difference 
in the slope of li is not obvious; however, the deviation in 
the ratio o f the le ngth to the width appears to be significant. 
IZou et ^ (120141 ') suggested using a render frame with the 
lighting of the radar model and combining multiple optical 
images using computer graphics methods, which may yield a 
likely explanation for the higher value of the length-to-width 
ratio obtained for the optical images acquired by C hang’e-2. 

_ Using the optical images from Chang’e-2, IBu et aP 

1 I 2 OI 5 II rotated the radar-derived model and retrieved an 


orientation with respect to the graphical frame described 
by direction cosine angles of (126.13 ± 0.29°, 122.98 ± 
0.21°, 126.63 ± 0.46°). In accordance with the dehnition of 
the cosine angles given bv iBu et al.l ll2015h . we calculate this 
set of cosine angles to be (130.3± 1.0°, 134.78±1.0°, 106.95± 
1.0°). The two sets of results are quite similar to each other. 
However, we should note that the results presented here ne¬ 
glect the attitude of the camera in the inertial frame. 

2.2.2 Euler angles 

Because of the fixed direction of the camera’s optical axis, 
Chang’e-2 maintained a nearly constant attitude throughout 
the shooting process. Figure 3 depicts the spacecraft’s body- 
fixed frame, where I is the direction of the camera’s optical 
axis and the corresponding unit vector in the spacecraft’s 
body-fixed frame is (—0.06976, 0.9976, 0). The location rela¬ 
tionship represents a transformation from Chang’e-2’s body- 
fixed frame to the graphical frame shown in Figure Ic. In 
combination with the attitude information of the spacecraft, 
we can establish a relationship between the graphical frame 
and the inertial coordinate system. As Figure Ic shows, the 
left portion of Toutatis is blocked by the solar panel on the 
spacecraft. The unit vectors of the x and y axes, which 
represent the horizontal and vertical axes in the graphi¬ 
cal frame, are determined to be (0.824,0.540,-0.173) and 
(0.344, —0.233, 0.91 0), respectively, in t he J2000 equatorial 
coordinate system dHuang et al.ll2013bl '). 

By merging the rotated radar-derived shape model and 
the optical image. Figure Id illustrates the similarities and 
differences between the two models. In combination with 
the conversion relationship among the asteroid’s body-fixed 
frame, the graphical frame, and the inertial coordinate sys¬ 
tem, we can obtain the matrix describing the transformation 
from the J2000 ecliptic coordinate system to the asteroid’s 
body-fixed system expressed in terms of 3-1-3 Euler angles 
as below: 

R = i?,(42.2°)i?427.6°)i?,(-20.1°)r, (1) 

where and Rz are the standard rotation matrices for 
right-handed rotations around the X and Z axes, respec¬ 
tively. The coordinate transformation and the corresponding 
Euler angles will be briefly introduced in the next section. 
The orientation of the principle axis is then obtained with re¬ 
spect to the attitude of the radar model, which is estimated 
to be (2 49.874:1°, 62 .43±1°) in the J2000 ecliptic coordinate 
system (IZhao et al.ll2014^ . Corresponding errors arise from 
the matching process, the uncertainties of the radar-derived 
shape model and the attitude uncertainties of the spacecraft 
(iHuang et al.ll20i^ ~). 

3 NUMERICAL MODELS 
3.1 Dynamical Model 

The rotation matrix shown in equation (1), which is com¬ 
posed of the three 3-1-3 Euler angles a = {a,P,^), maps 
the information for conversion from the inertial coordinate 
system to the body-fixed frame, as shown in Figure 4. The 
body-fixed frame is generated by applying the following ro¬ 
tation sequence of the yaw, pitch and roll angles: a rotation 
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Figure 2. Properties of an optical image of Toutatis. a: The red line and green line represent the lengths of the image of the asteroid 
along the li and I 2 directions, respectively; the ratio between these lengths is one of the properties that is used to characterize the image. 
The blue line represents the horizontal line, which can be used to determine the slope of li. b and c: The area within the green circle 
represents an obvious characteristic topological feature that can be used to derive information concerning the rotation about li. 



Figure 3. a. Graphical frame, b. Spacecraft’s body-fixed frame. 


of a around the z axis, then a rotation of /3 around the x 
axis, and finally, a rotation of 7 around the z axis from the 
inertial frame. 

The Euler angles describe the orientation of a rigid body 
in space at a specific time, and their variations represent spin 
states. Let the vector uj define the instantaneous rotational 
velocity in the body-fixed frame. Then, the set of kinematic 
differential equations for the 3-1-3 Euler angles is as follows: 


1 

sin/3 


sin 7 

cos 7 sin /3 
— sin 7 cos P 


cos 7 0 \ 

— sin 7 sin /3 0 1 

- cos 7 cos P sin P ) 

= [B{a)]a 


( 2 ) 


The time derivative of the Euler angles encounters a 
singularity at either /3 = 0° or /3 = 180°, which may cause 
a and 7 to rotate in the same plane. 

As derived from the Euler equation, Euler’s rotational 
equation of motion describes the time derivative of the an¬ 
gular velocities: 

[/]a3 = — [ii)][7]a;-I-L , (3) 


where the moment-of-inertia matrix [7] is constant, symmet¬ 
ric and given in the body-fixed frame. It has dimensions of 
[3 X 3] and can be defined in terms of six quantities. L repre¬ 
sents the external torques acting on the dynamical system, 
and the matrix [cj] has the following form: 


[<i] = 


0 — cJa 0^2 

CJ3 0 — uj\ 

—U2 UJl 0 


( 4 ) 


Then, the time derivative of the angular velocities is 
expressed as follows: 

^ = [i]-\-[m^+L) ■ ( 5 ) 


To calculate the external torque exerted by a spherical 
body, we assume that r* is the vector of an infinitesimal 
mass element (dm) of the asteroid relative to the origin. 
Figure 5 presents a schematic diagram of the effect exerted 
by a spherical perturbing body on an irregularly shaped rigid 
body. According to the definition of angular momentum H, 
we have 


_ _d 
dt dt 


dm r X f' = 


r X (dFc/dm — B)dm 


( 6 ) 


where B is the acceleration of the origin in the inertial co¬ 
ordinate system and dFc is the gravitational attraction ex¬ 
perienced by an arbitrary mass element. When the center of 
mass of the rigid body is chosen as the origin O, we have 
f dmr X B = 0. Therefore, th e torques exerted o n the rigid 
body are written in the form JSchaub et ahlfioogl l 


Lo 



r X dFc 


( 7 ) 
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Equator of 
Toutatis 


Figure 4. Coordinate system transformation relationship in terms of 3-1-3 Euler angles. Subscripts i and b indicate the inertial coordinate 
system and the body-fixed frame, respectively. The three angles a, and 7 form a 3-1-3 set of Euler angles. 


Otherwise, if the vector of (dm) relative to the center of 
mass is f),, equation (6) has the following form: 


dH /* -» 

= I r X (dFc/dm — B)dm 


= / (f7 - r^m) X (dFc/dm - B)dm , 


( 8 ) 


where Tcm = (xcm, Vcm, Zcm)'^ is the vector of the center of 
figure with respect to the center of mass (the COM-COF 
offset). Thus, equation (6) can be expressed as follows: 

L= r X dFa + Mvcm x B , (9) 

J M 

where M is the total mass of the asteroid. 

Let r be the position vector of the spherical body rel¬ 
ative to the asteroid. R = ^ — r denotes the position of 
a body element measured from the spherical body. Conse¬ 
quently, dFc has the following form: 

, - GMs g . 

dFa =- —^Rdm , 

Rr 

where Mg is the mass of the external spherical body. Hence, 
the integral in equation (9) can be expressed as 


( 10 ) 


( 11 ) 


L = -GMg r X ^dm -f Mf^m x B . 

Jm R 

As Figure 5 shows, under the assumption that r ^ r', 
the hrst-order expansion of has the following form: 


^3 1 ^ ^2 J ^3 ^ ^3 


( 12 ) 


Substituting (12) into equation (9) allows the moment 


from the perturbing body to be divided into two compo¬ 
nents: 


L — Li -t- Z/2 , 


( 13 ) 


where 


Li 



■^dm -I- MFcm X B , 


L2 = 


3GMs 



r){r X r)dm . 


(14) 

(15) 


If the origin O is the center of mass of the irregularly shaped 
rigid body, then we obviously have rym = 0 and Li = 0. 
Otherwise, truncated at the lowest order, we have 



(16) 


If perturbing bodies other than the Sun are considered, 
then equation (16) should incorporate their contributions. 
Substituting equation (16) into equation (14), we obtain 


Li = GMa {r + Vcm) X ^dm 




(17) 


Essentially, the torques imposed by the COM-COF off¬ 
set are far too small to contribute in the first-order approxi¬ 
mation. Furthermore, the change in angular momentum that 
results from Li is not an intrinsic feature of the rigid body; 
rather, it is determined solely by the selection of the base 
poi nt when solving for the rotational dynamics. 

iTakahashi et al.l (l2013f ) expressed the second-degree po¬ 
tential as follows: 
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perturbing body 



Figure 5. External torque of a spherical body on an irregularly shaped rigid body. 




(18) 


where It is the trace of [/]. The partial derivative of the 
potential can be derived as follows: 


dUi 

dr 


3G ^ 



(19) 


therefore, the deduced moment can be expressed in the fol¬ 
lowing form: 


L 2 


3GMs 


[r][I]r . 


( 20 ) 


Our second-degree expansion of the moment given in equa¬ 
tion (15) is consistent with the results represented by equa¬ 
tion ( 20 ), which indicates that tidal torques caused by the 
oblateness of the rigid body. 


3.2 Numerical Method 


The least-squares and multiple shooting methods are used 
to fit the observational data and to simulate the propagation 
of the rotational parameters. The dynamical equation of the 
state vector, which is composed of three Euler angles and 
three angular velocities, is written as follows: 

X = F{X,t) , (21) 

where X = (a,/3, 7 , cua,, cuz)^; a, jS, and 7 are the 3-1-3 
Euler angles, and uji, uj 2 , and ljs are the three components of 
angular velocity. The previous results for Iw, Izz, Ixv, 
Ixz, and lyz are used in our calculations (iTakahashi et ah] 
I 2 OI 3 II . F is the function that represents the time derivative 
of X and the dynamical model. In the first-order approxi¬ 
mation, the dynamical model has the following form: 

X = A{X,t)X{t) , (22) 


where A{X,t) is the dynamical matrix. The transition ma¬ 
trix > 1 > is defined as follows: 


<&(tl,t2) 


dx{ti) 

dx{t2) 


(23) 


Based on the nominal orbit, we establish the relation¬ 
ship between the state vectors at two specific times, and 
t 2 , as follows: 


X{t2) = <l){ti,t2)X{ti) . (24) 

Then, the transition matrix can be obtained by inte¬ 
grating the derivate equation: 

= A{t,X)(j){to,t) . (25) 

The observational equation for the measurements Y is 

Yi = Z{Xi,ti) + ei , (26) 

where Z represents the observational model, the subscript 
i indicates the sequence of the observational data and e 
represents the observational uncertainties. In the first-order 
approximation, the partial derivatives of the observational 
quantities with respect to the variables form the relation¬ 
ships between them, and we have 

dZ dZ 

Yi = \x—7:iXi -j- ei = ti)X{to) Ci , (27) 

where ti is the observation epoch of the data set. The 
least-squares method adopted herein is similar to that of 
ITakahashi et al.l l| 201 ij l. The cost function is defined as fol¬ 
lows: 

f)7 

+i(X(to)-X(fo))^p-' 

(X(to)-X(fo)) , 

(28) 

where W is the weighting matrix, and P is the covariance 
matrix. The bars over X and P represent the a priori values 
deduced through estimation or from previous results. The 
modified differential equation that computes the correction 
to the variables can be expressed as follows: 

AX{to) = 

C^H^WH + p-^y^'^H^WY + p-^X{to)y^ . (29) 

The RKF78 integrator can be used to integrate the rota¬ 
tional equations. The initial step size is set to approximately 
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10-9 o Euler angles and 10“® ° jday for the rotational 

velocities. In the integration, an adaptive step size is used 
to numerically solve the equations. The truncation errors for 
the Euler angles and rotational velocities are set to 10“® ° 
and 10~^ °/day, respectively. 


4 RESULTS 

To better understand the rotational dynamics of Toutatis, 
we performed a large number of numerical simulations based 
on Chang’e-2’s observations and ground-based radar mea¬ 
surements, using our dynamical models described above. In 
the following, we will present the major outcomes for the 
spin states, rotational period and variation in angular mo¬ 
mentum of Toutatis. 

In this work, the initial variables for Toutatis’ spin state 
that were used in ou r num erical simulations were adopted 
from iTakahashi et al.l (l2013l ~) for the epoch to (17:49:47 UTC 
on 9 Nov 1992). The orientation of the angular momentum 
and the rotational periods were also calculated in the study. 
In our calculation, we balanced the weights of the optical 
data based on the uncertainties of the observations. Conse¬ 
quently, we derived innovative solutions for the spin state of 
Toutatis, which are summarized in Table 3. 

The residuals of the Euler angles (34 sets) and angular 
velocity (33 sets) were normalized with respect to the maxi¬ 
mum radar observational errors (15° for the Euler angles and 
10°/day for the angular velocities) and are shown in Figure 
6. Because the previous prediction of the orientation at the 
Chang’e-2 flyby epoch based on the radar-derived results dif¬ 
fers from that observed by Chang’e-2, the use of the optical 
data might have degraded the convergence of the simulation 
algorithm. Therefore, the magnitude of the residuals is larger 
than that observed in the results from the radar data (Taka- 
hashi et al, 2013). The residual errors in the simulations were 
normalized with respect to the observational uncertainties. 
Because of the inconsistency of the observational data, the 
magnitudes of the residuals are slightly higher than those of 
the previous results. However, all deviations lie within the 
3a region. The largest bias is found in the roll angle, which 
exhibits a remarkable difference between the prediction ob¬ 
tained from the radar measurements and the authentic spin 
state of Toutatis that is directly indicated by Chang’e-2’s 
observations at the flyby epoch. 


4.1 Spin States 

According to the numerical re sults derived from radar obser¬ 
vations collected before 2008 dBusch et al.ll201^ ). we consid¬ 
ered a render effect for Toutatis and generated a predicted 
imaging outcome prior to Chang’e-2’s flyby, as shown in Fig¬ 
ure 7a llBusch et ahlboi^ : IZhao et al.ll20143 l. In addition, 
based on the images acquired by Chang’e-2, we corrected 
the attitude of Toutatis by rotating the radar-derived shape 
model (see Section 2.2.1) to search for a good match with 
the Chang’e-2 images acquired at the flyby epoch (Fig. 7b), 
which provide the only space-borne optical data regarding 
Toutatis’ orientation. Furthermore, the present simulations 
yielded another solution for Toutatis’ attitude during the 
near-Earth flyby in 2012. Figure 7c shows the outcomes de¬ 


rived from our rotational model using space- and ground- 
based observations. 

In comparison with the results obtained from the op¬ 
tical images (Fig. 7b), the radar-derived results (Fig. 7a) 
exhibit a dramatic deviation in the roll angle; hence, these 
results yield a different profile of the asteroid. The simula¬ 
tion results derived from our dynamical model (Fig. 7c) dif¬ 
fer from those of Figure 7b with a pitch angle bias of within 
20°. Thus, we may safely conclude that our outcomes repre¬ 
sent a good improvement in the understanding of Toutatis’ 
spin state. 

The orientation of the long axis in the inertial frame 
likely reflects the precession of Toutatis. Based on the dy¬ 
namical model of rotation, we calculated the variation in the 
direction of the long axis. Figure 8 shows the trajectories of 
the long axis with respect to the J2000 ecliptic coordinate 
system in a unit sphere over the past two decades. The mo¬ 
tions of the long axis are projected onto the X-Y, X-Z and 
Y-Z planes (see Figs. 8a, 8b and 8c, respectively). The fig¬ 
ure reveals that the long-axis motion of the asteroid has 
remained ellipsoidal in the X-Y and Y-Z planes, whereas it 
has rectilinearly precessed in the X-Z plane. All curves lie 
outside the ecliptic plane, implying that the small lobe of 
Toutatis is always located above the large end from a view¬ 
point close to the ecliptic. 

Moreover, the orientation of the center axis of preces¬ 
sion of Toutatis can be approximately determined from Fig¬ 
ure 8, and the derived spherical coordinates can be estimated 
to be (—0.2°, 54.6°) in the ecliptic coordinate system. As 
a result of Toutatis’ clockwise rotation and precession, the 
center axis of precession points nearly along the opposite 
direction to the angular momentum (see Section 4.3). The 
amplitude of the precession is approximately 60°, which may 
shed light on the significantly different attitudes of the as¬ 
teroid that have been observed from Earth. 

4.2 Angular Momentum 

Considering the radar-derived shape mod el and the com¬ 
ponen ts of the inertial matrix inferred by iTakahashi et al.l 
ll2013tl . we will now explore the angular momentum of Tou¬ 
tatis induced by various external gravitational torques. Fig¬ 
ure 9 shows the variations in the external gravitational 
torques acting on the spin state of Toutatis from 1992 to 
2012. The solar torque is on the order of 10“^° — 10“^, indi¬ 
cating that its value is 2-3 orders higher at the perigee than 
at the apogee. Its periodic variation is clearly associated 
with Toutatis’ orbital period. The variation tendencies of the 
gravitational torques arising from the Earth and Moon are 
similar, as shown by the red and blue curves, respectively. 
There is a 10“^ difference in the orders of these torques be¬ 
cause of the magnitudes of the masses of the bodies from 
which they arise. The periods of both torques are consistent 
with that of the black curve because of Toutatis’ resonance 
orbit with the Earth. At present, Toutatis is also in a 3:1 
mean motion resonance orbit with Jupiter dWhipple et al.l 
I 1993 II : thus, one entire period of the torque induced b y 
Jupiter is displayed by the green curve (iBusch et al.ll20f3 l. 

Based on the rotational dynamical equation and the in¬ 
tegrated orbit, the overall influence of these external torques 
on the variation in the magnitude of Toutatis’ rotational 
angular momentum from 1992 to 2012 was normalized with 
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Table 3. Spin-state parameters of Toutatis derived from our numerical simulations. 


Property 


Value 


Simulated solutions at to 

Results at flyby epoch 
Orientation of angular momentum 
Rotational/precession period 


o; = 147.5°,^ = 63.9°,7 = 241.5° 
uji = 14.5°/day, aJ 2 = 33.7°/day= 
a = -3.65°,/3 = 43.62°,7 = 24.7° 

\h = 180.2+°;2°,/3^ = -54.757;}®° 
5.38 days, 7.40 days 


—9S.6°/day 



Figure 6. Residuals of rotational parameters with respect to the maximum radar observational errors ( 15° for the Euler angles and 
10°/day for the angular velocities). Left panels: residuals of the Euler angles. Right panels: residuals of the angular velocity components. 



Figure 7. Comparison of Toutatis’ orientation, a: The left panel shows the results derived from iTakahashi et al.l l|2013h ’s work, with an 
uncertainty of dozens of degrees, b: The middle panel shows the results calculated from optical images acquired by Chang’e-2 and by 
rotating the radar shape model, c: The right panel presents the outcomes derived from our rotational model using space- and ground-based 
observations, which are close to those obtained from the optical images acquired by Chang’e-2. 
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Figure 8 . Trajectories in the J2000 ecliptic coordinate system of the long axis of Toutatis in a unit sphere. Panels a, b, and c show the 
trajectories in the X-Y, X-Z, and Y-Z planes, respectively, and panel d shows the motion of the long axis in 3D. 



Time (years) 


Figure 9. Variations in external gravitational torques. The black, red, blue, and green curves represent the torques arising from the 
Sun, the Earth, the Moon, and Jupiter, respectively. 


respect to the initial magnitude Ho llTakahashi et al.]l2013ll . 
as shown in Figure 10a. The terrestrial tidal torque (see the 
red curve in Fig. 10b) causes a considerable change in angu¬ 
lar momentum when the asteroid approaches Earth at the 
perihelion or during the Earth flyby that occurs every four 
years. The most significant change, with a variation in angu¬ 
lar momentum magnitude on the order of 0.03%, occurred 
in 2004 as a result of Toutatis passing the Earth within 4.02 
lunar distances. Similarly, the tendency of the effect of the 
lunar torque is consistent with that of the terrestrial torque, 
as shown in Figure 10c. The solar tides always have a pre¬ 


dominant influence on the rotational variation. However, the 
terrestrial tides also play an important role in the variation 
in angular momentum during Toutatis’ regular nearby vis¬ 
its to Earth. Figure lOd shows the influence on the angular 
momentum exerted by Jupiter. The order of magnitude of 
this effect slightly changed after the 2004 near-Earth flyby, 
and the amplitude continually remains lower than that of 
the terrestrial torque. 

As our simulation results indicate, the angular momen¬ 
tum orientation of Toutatis is determined to be described 
by (Air = 180 . 2 |J]g 3 o and Ph = —54.75^q'5^qo) and has re- 
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Figure 10. Variation in the angular momentum of Toutatis from 1992 to 2012. 


mained nearly unchanged in space over the past two decades. 
Figure 11 shows the variations in Toutatis’ angular momen¬ 
tum orientation from 1992 to 2012 in the J2000 ecliptic 
frame. The amplitude of this change is shown to be less than 
one degree in both longitude and latitude. Jumps in the an¬ 
gular momentum orientation occur at the perihelion of each 
orbit. A small change in behavior is evident in Figure 11b as 
a result of the 2004 near-Earth flyby, consistent with Figure 
10b. The variation with solar distance that is apparent in 
Figure 11c and lid indicates that the solar and terrestrial 
torques predominantly affect the rotational motion of the 
asteroid. The misalignment of the curves in Figure lid is 
also a result of the 2004 near-Earth flyby. Figure He and 
Ilf show the angular momentum orientations of 33 sets of 
radar observations. Compared with the numerical results, 
the observation data fall within a reasonable error range, 
with the exception of a few points at large bias. 

4.3 Rotational Period 

The motions of the short or middle axis reflect the status 
of the rotation about the long axis, whereas the long axis’ 
motion represents precession. To calculate the two periods 
associated with the spin states of Toutatis, we determined 
the latitudinal variations of the asteroid’s long and middle 
axes in the J2000 ecliptic frame, as shown in Figure 12. We 
applied Fourier transform to analyze the periods of the two 
oscillation parameters and found that they are 5.38 days for 
the rotation about the principal axis and 7.40 days for the 
precession of the principal axis. Th ese results are in good 
agreement with the previous results dOstro et al.ll 19991 ). 

Let Px and Pz indicate the latitudes of the asteroid’s 
long and middle axes, respectively, in the J2000 ecliptic coor¬ 
dinate system. Figures 12 and 13 show the latitudinal varia¬ 
tions of these axes during the 1992 and 1996 flybys, respec¬ 


tively. Additionally, the numerical results (represented by 
dotted lines in Figs. 12 and 13) that were calculated from 
our dynamical model are found to be in good agreement 
with the radar observations wit hin the error bars (marked 
by stars) llTakahashi et al.]l2013l ). as listed in Table 4. Hence, 
we may conclude that the orientation parameters of Tou¬ 
tatis obtained from our investigation are very reliable. This 
evidence provides further confirmation that our proposed 
rotational model can be used to correctly evaluate the spin 
status of Toutatis or other asteroids. 


5 CONCLUSIONS AND DISCUSSIONS 

In this work, we apply the observations collected during 
Chang’e-2’s outbound flyby to model the rotational dynam¬ 
ics and determine the spin state of Toutatis. Based on flyby 
images, we utilize the radar-derived shape model to calcu¬ 
late Toutatis’ orientation at the flyby epoch. In addition, we 
estimate the 3-1-3 Euler angles to be —20.1°±1°, 27.6°±1°, 
and 42.2° ± 1°, respectively. Consequently, our results have 
greatly improved the estimation of the orientational param¬ 
eters of Toutatis with respect to the previous predictions. 

In combination with ground-based observations, we in¬ 
vestigated the evolution of the spin parameters using nu¬ 
merical simulations. In addition to the solar and terres¬ 
trial torques, the tidal effects arising from the Moon and 
Jupiter are extensively considered in our dynamical model. 
The magnitude and influence of these gravitational torques 
were analyzed in this work. The solar tide appears to always 
be the dominant torque acting on the angular momentum 
of Toutatis. Furthermore, the contribution to the external 
gravitation torque due to the COM-COF offset appears to 
be negligible in the first-order approximation. We also found 
that the closest near-Earth flyby, at 4.02 lunar distances, re- 
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Figure 11. Variations in Toutatis’ angular momentum orientation from 1992 to 2012. Panels a and b show the variations in longitude 
and latitude versus time. Panels c and d show the change in longitude and latitude with the distance of Toutatis from the Sun. Panels e 
and f show the corresponding results from radar observations. 


suited in a 0.03% change in the magnitude of the angular 
momentum of Toutatis. The dynamical influence exerted by 
Saturn on the angular momentum was also assessed in fur¬ 
ther simulations and found to be approximately 10 ^ lower 
than that of Jupiter. Hence, we can safely conclude that Sat¬ 
urn plays a less important role in the variation of Toutatis’ 
angular momentum. 

The attitude at the Chang’e-2 flyby epoch that was de¬ 
rived from the numerical simulations yielded a better ap¬ 
proximation to the optical results than that previously ob¬ 
tained from radar data alone. The largest deviation in the 
Euler angles is observed in the pitch angle, with a bias of 
less than 20°. The uncertainties corresponding to observa¬ 
tional uncertainties and data processing error were consid¬ 
ered in the simulations. The inconsistency in the different 
types of observational data may have led to higher residu¬ 
als compared with previous results. Simulations based solely 
on radar observations were also performed, and the corre¬ 
sponding rms magnitude was much lower. However, the re¬ 
sults obtained using a higher-accuracy dynamical model and 
a combination of the various types of observations yielded a 
good result that is highly consistent with the optical images 
acquired during the Chang’e-2 flyby in 2012. 

The precession of Toutatis was investigated by consid¬ 
ering the motion of its long axis. The behavior of the orbit 


in the inertial frame was found to be circular, with a center 
axis pointing along (—0.2°, 54.6°). The precession amplitude 
was estimated to be up to 60°, which may be responsible for 
the significantly different attitude of the asteroid as observed 
by ground-based facilities. Moreover, by exploring the mo¬ 
tions of the long axis and middle axis, we determined the 
rotation period of Toutatis using Fourier analysis. The two 
major periods were found to be 5.38 days for the principal 
axis rotation and 7.40 days for the prece s sion, i n agreement 
with the results reported bv lOstro et al.l lll999ll . 


Toutatis’ angular momentum orientation was deter¬ 
mined to be described by Xh ~ 180 . 2 ^q’ 3 o and Ph = 
—54.75^q jqo , indicating that it has remained nearly un¬ 
changed for the last two decades. Because of the increasing 
magnitudes of the solar and terrestrial torques, tiny jumps 
in the angular momentum orientation occur at perihelion in 
each orbital period. However, the dynamical effects caused 
by the near-Earth flyby in 2004 slightly changed the lat¬ 
itude of Toutatis’ angular momentum orientation. Hence, 
our simulation results are in good agreement with previous 
radar observations. In a word, based on the combination of 
Chang’e-2’s observations and radar data, our investigation 
offers an improved understanding of the rotational dynamics 
of Toutatis. 
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Figure 12. Latitudinal variations of Toutatis’ long axis Pz and middle axis in the J2000 inertial coordinate system in 1992. 
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Figure 13. Latitudinal variations of Toutatis’ long axis j^z and middle axis in the J2000 inertial coordinate system during the 
near-Earth flyby in 1996. 
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